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We examine some aspects of the continuous photodetection model for photocounting processes in 
cavities. First, we work out a microscopic model that describes the field-detector interaction and 
deduce a general expression for the Quantum Jump Superoperator (QJS), that shapes the detector's 
post-action on the field upon a detection. We show that in particular cases our model recovers the 
QJSs previously proposed ad hoc in the literature and point out that by adjusting the detector 
parameters one can engineer QJSs. Then we set up schemes for experimental verification of the 
model. By taking into account the ubiquitous non-idealities, we show that by measuring the lower 
photocounts moments and the mean waiting time one can check which QJS better describes the 
photocounting phenomenon. 
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I. INTRODUCTION 
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The subject of quantum measurements is as old as the very foundation of quantum mechanics. For a long time 
the scheme of sudden state reduction, proposed by von Neumann, has been prevalent. He conjectured that the 
measurement of an observable on a system entails its state reduction to one of its eigenstates, or shortly, a sudden 
change of the system state by projection. However, on probing an electromagnetic (EM) field state through a 
photocount process, the photons are detected and counted one by one, a photon entering a photomultiplier tube 
provokes a burst of electrons (a photocurrent) which is viewed as originating from that single photon. It is then 
registered and counted. A sequence of bursts in a given time interval is associated to the photocount process. So, 
the determination of the field state is not achieved by an instantaneous projective measurement, but it takes some 
time t to count a sequence of photons, whose statistics gives information about the field state. A classical theory 
describing this process was proposed by Mandel [l[ . Quantum photocount theories were developed by several authors 
@, H II H Hi (see the review [7J for more references). These theories rely on the assumption of an instantaneous 
£7^" measurement of m photons, independently of the duration of the detection. However, actually, photons are counted 
sequentially, one by one, and the time intervals between counts is irregular and uncontrollable. 

For describing more realistically a sequential photocount events in an ideal closed cavity, Srinivas and Davies 
' (SD) Q proposed an approach based on the concept of continuous measurement. Their scheme allows calculating 
various statistical functions, that can be compared with experimental outcomes, such as the probability for counting 
any number k of photons in a time interval t and different coincidence probability densities. The SD model takes into 
account a "back action" of the photodetector on the state of field and gives the conditioned field state, i.e., the field 
state just after a given sequence of photocounting events. A progress in understanding the physical meaning of the 
axiomatic SD model was achieved due to studies [3, El, E2> EH (for other references see m 113). 

Continuous photodetection model (CPM) is extensively discussed in the literature fA foL llll Il5l ITfij. so we shall 
mention only its main properties. The model, also referred as a theory, describes the field state evolution during 
the photodetection process in a closed cavity and is formulated in terms of two fundamental operations, assumed to 
represent the only events occurring at each infinitesimal time interval. (1) The one-count operation, represented by 
the Quantum Jump Superoperator (QJS), describes the detector's back-action on the field upon a single count, and 
the trace calculation over the QJS gives the probability per unit time for occurrence of a detection. (2) The no-count 
operation describes the field non-unitary evolution in absence of counts. 
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If one sets the formal expressions for these operations, all possible outcomes of aphotocounting experiment can be 
predicted. For instance, the photocounts [1, [nj, EH and the waiting time (WT) [l7], [H, [l9|, [2(| statistics are among 
the most common quantities studied both theoretically and experimentally. Moreover, CPM conferred a new step 
in photodetection theories by allowing one to determine the field state after an arbitrary sequence of measurements, 
thus creating the possibility of controlling the field properties in real time experiments |26l [2a . |53| . 

Actually, the QJS is the main formal ingredient within the theory, since it also dictates the form of the no-count 
superoperator Q. Two different models for the QJS were proposed ad hoc. The first was proposed by Srinivas and 
Davies (SD) in the original paper Q (we call it SD-model) as 

JsDP = Tlapa\ (1) 

where p is the field statistical operator, a and a' are the usual bosonic annihilation/creation operators and TZ is 
roughly the detector's ideal counting rate [1, Hlj]. From the very beginning the authors [|[ denounced the presence of 
some inconsistences when Js d is employed for describing a real photodetection process, this point was also appointed 
in [HI. Nevertheless, this QJS is widely used in the literature p], 0, H H Hf H M, E3, IH, [U, HIM ill. 

The other proposal 0, [H| assumes for the QJS an expression written in terms of the ladder operators E- = 
(a) a + l) 1 2 a and E+ = E\_ (also known as exponential phase operators [34L [35l l36l l37l. l38l| ) 

J E p = 1lE-pE + . (2) 

In [2l| we called E-model such a choice, to distinguish it from the SD QJS |T]). We note that TZ may be different for SD- 
and E- models, but the above notation will not cause confusion in this paper. Besides eliminating the inconsistencies 
within the SD-model, the use of the E-model leads to different qualitative and quantitative predictions for several 
observable quantities. 

In section |TI] we present a microscopic model for the detector assumed to be composed of a sensor (2-level quantum 
object) and an amplification mechanism (macroscopic thermal reservoir). In section [IIII we compare our model's 
predictions concerning photodetector properties with experimental data and show that the QJSs ((T|) and ^ are 
particular cases of a general time-dependent transition superoperator, each one occurring in a particular regime of the 
detector experimental parameters [39l. |40||. Moreover, we point out that by manipulating detector's parameters one 
could engineer the form of the QJS, thus changing the dynamics of the photodetection, as well as the field state after 
a sequence of measurements. 

A way to check the validity of CPM and to decide which QJS better describes the phenomenon in practice can 
be accomplished through photodetection experiments in a high finesse cavity by comparing experimental outcomes 
with theoretical predictions. However, real detectors and cavities are far from ideal. So in section HVl we include the 
main non-idealities [quantum efficiency (QE) and dark counts] into the CPM and deduce general expressions for the 
photocounts and the WT distributions. As a practical application, in section [V] we give some experimental hints to 
decide which QJS model actually prevails in a photodetection experiment. Section [VII contains a summary and the 
conclusions. 



II. MICROSCOPIC MODEL OF 2-LEVEL PHOTODETECTOR 



We model the photodetector as constituted of two parts: the sensor and the amplification mechanism (AM). The 
sensor is a two-level quantum object (atom-like) with resonant frequency luq, interacting with the monomodal EM 
field of frequency lo. It has the ground \g) and the excited e) states, so we describe it by the usual Jaynes-Cummings 
Hamiltonian [41| 

H = -luqcto + + gaa+ + g*aV_, (3) 

where g (assumed to be real, since only its absolute value enters the final expressions) is the sensor-field coupling 
constant, and the sensor operators are o~q = |e)(e| — |<7)(g|, c + = \e)(g\ and tr_ = |ff)(e|. The interaction described by 
Hamiltonian ((3j allows a coherent exchange of excitations between the sensor and the field - the Rabi oscillations. 

Upon absorbing a photon the sensor initially in \g) makes a transition to |e), and after some time it decays back, 
emitting a photoelectron into the AM; after that, the detector is ready for the next photodetection. The AM is a 
complex macroscopic structure (e.g. photomultiplicr tube) that somehow amplifies the photoelectron and originates 
some observable macroscopic effect, giving rise to the click of the detector. In order to describe general features of 
the AM independent of the type of the single photon detector (SPD), we model it as a macroscopic thermal reservoir 
with a mean number of intrinsic excitations n (due to the effects of temperature and internal defects). Thus, the 



3 



whole system field-SPD unconditioned time evolution (when the detector is not monitored [14|]) is described by the 
master equation [401 ] 

p T = - [Ho,p T ] - jn (<j-<T + p T - 2<T + p T <7- + p T <j-<7 + ) 
i 

—7 (n + 1) (cr + er_p T - 2cr_p T cr + + p T a + a _) , (4) 

where 7 is the sensor-AM coupling constant. 

According to CPM, the trace of QJS applied on the field density operator gives the probability density p(t) for the 
photodetection, i.e. emission of a photoelectron at time t, given that at time t = the detector-field system was in 
the state 

Po = \9)(g\ ®P, (5) 

where p is the field statistical operator. Microscopically this means that initially the detector is in the ground state; 
then, during the time interval (0, t) the sensor interacted with the field and it could have absorbed a photon, doing 
a transition \g) — > |e). So p(t)A£ is the probability of the sensor decaying back to \g) during the time interval 
(t, t + At) and simultaneously emitting a photoelectron that will lately originate one click. Here, the emission of the 
photoelectron is our interpretation of how the detector operates, and this phenomenon does not appear explicitly in 
the formalism. 

Following the quantum trajectories approach (T^Ji p(t) is calculated as 

p(t) = Tr F _ D [MtPo] , (6) 

where Utp represents the evolution of the field-SPD system from initial state p at time t = to the time t without 
detections, and RUtPo stays for a click (instantaneous decay of the sensor) at the time t (in the trace, F stands for 
the field and D - for the detector). The sensor instantaneous decay |e) — > \g) is represented by the superoperator 

RPo = 27 (™+ l)o--Po cr +> ( 7 ) 

whose trace gives the probability density of such an event. In Eq. (7]) 7 stands for the sensor zero-temperature decay 
rate, and we included the term n + 1 because it is natural to assume that the rate of decays is proportional to the 
effective temperature of the detector (proportional to n). The complementary no-decay superoperator U t p — pjj(t) 
describes the non-unitary evolution of the field-SPD system during time interval (0, t) without clicks; it is the solution 
to the master equation ((4]) without the decay term ||7J): 

= -i(H ePu - Pu Ht) + 2 1 n<7 +Pu a_ 1 (8) 

where 

(ujo — vy) „ + ._ 1. . . 

H e = erg + uju + gacr + + ga ] o _ — ij{n + —). (9) 

Taking the partial trace over the detector variables in ([6]) one obtains the superoperator 

Z(t)p = Tr D [RU tPo \ , (10) 

which describes the back-action of the detector on the field upon one click - it is the transition superoperator [39[ 
and, as will be seen below, its time average defines the QJS. Moreover, the probability density for a count is simply 
p(t)=Tr[E(t)p}. 

In order to solve Eq. © we first do the transformation 

Pu = X t ~p v x\, X t = exp(-iH e t) (11) 

to obtain a simple equation for p v 

= 27rT(T + p [/ CT_, <7_|_ = A_ t cr + A 4 , = a\, (12) 
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whose formal solution is 

Pu(t) = Po + 27" / di'a+(t')MOMO- (13) 



If one iterates Eq. ([13]) and substitutes the result into Eq. ([10]) one gets the following transition superoperator 

oo 

E(t)p= 2bg (l + n)J2 (^n) 1 (t)p, (14) 
z=o 



where 



6 = 7/5, So(t)p = Mj> e Q (t) = (e\X t \g). (15) 



For I > 

E t (t)p= [ dtf- [ dtihpe], 0/(Mi,"- = (e\X t a + (t 1 )a + {t 2 )---a + (t l )\g). (16) 
Jo Jo 

As shown in [4p| it is enough to evaluate just the three initial terms in the sum (|14p . whose constituents are found 
to be 

O = -ie-^+W-Mn+l/2)tg A+i ^ g i = e - 7 t(W+l/2)- iw ftt Xft+i(i _ fc^-t^ 

2 = - ie -i*Cs+i/2)-*- (ft-i/2)t x . +i(t _ h)Sfi (tl _ taj^f-tajat, 



where 



C fi (t) = cos , = sin ( 7 tB fl /6) /B fl , Xfl = e~" t/2 [C ft (f) - , (17) 



= Vn + S 2 , 5={q- ib)/2, q=(uj - w) /$. (18) 

Substituting these expressions into Eq. l[T4|) . the transition superoperator turns out to be time-dependent, contrary 
to the standard definition of the QJS. So we evaluate the QJS as the time average of over the time interval T 
(to be determined later) during which the photoelectron is emitted with high probability 

Jp=^J o dtZ(t)p. (19) 

Considering the weak coupling (u),u>o 3> 7,9) for which the Jaynes-Cummings Hamiltonian |3|) and the master 
equation ^ are valid, and expressing the field density operator in Fock basis as 

00 

P= Pmn\ m )( n l (20) 

m,n— 

after the averaging in ([19]) the off-diagonal elements of Jp vanish due to rapid oscillations of the terms exp(±«w£). 
Therefore, we are left only with the diagonal elements of J p. Applying the superoperators on the density matrix 
as in ([T4"|) and evaluating Eq. (TT5]) we obtain 

00 

■b = Y.Pnn [nJ { n B) \n-l)(n-\\ +j( D )\n)(n\ + (n + 1) J^ } |n + l)(n + 1| + • • •] , (21) 

n=0 

where the explicit expressions for the n-dependent functions Jr t B \ Jn and are given in [io| . 

The QJS pTjl contains an infinite number of terms, so after a click the initial field state p reduces to a mixture of 

( B) 

different states, each one with its respective probability. The first term, with coefficient Jn , takes out a photon from 
the field, so it represents a click preceded by a photoabsorption - we call this event a "bright count" . The second term, 
dependent on J„ (proportional to n, quite small as will be shown below), does not subtract photons from the field 
but only modifies the relative weight of the field state components - it represents a "dark count" , when the detector 
emits a click due to the amplification of its intrinsic excitations. All further terms in Eq. (|2ip are proportional to n l , 

I > 2; they describe emissions of several photons into the field upon a click, so we call the first of these term, J„ , 
the "emission term" . There are many different phenomena that give rise to dark counts [i^, IH, S3] , our model takes 
into account only those causing the sensor's ground-excited state transition. Since the sensor's state depends on its 
interaction with the field, the dark counts modify indirectly the relative weight between the field state components - 
that is why they depend on n, what is not obvious at first glance. 
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FIG. 1: Signal-to-Noise ratio as function of b at resonance FIG. 2: Bright counts rate as function of wavelength of the 
(A = 500 nm) and far from it (A = 1500 nm) in the inset, field for different values of b, showing that TZb increases pro- 
The estimated breakdown value is &s ~ 380. portionally to b. The resonant wavelength is Ao = 500 nm. 



III. COMPARISON TO EXPERIMENTAL DATA AND QJS ENGINEERING 

Now we compare our results regarding SPD properties with available experimental data. Experimentally [4r3 | the 
dependencies of bright and dark counts rates are set as functions of the light wavelength and the detector's "bias 
parameter" (BP). In the BP we englobe such quantities as bias voltage, bias current and other physical quantities 
the experimenter adjusts in order to achieve simultaneously the highest Signal-to-Noise ratio S and bright counts 
rate. S = TZb/TZd is the ratio between the bright, TZb, and dark, TZd, counts rates. When one increases the BP, the 
bright counts rate increases while S remains unchanged until the breakdown value of BP, when S starts to fall rapidly 
as function of the BP. So most detectors usually operate near the BP breakdown in order to achieve the optimal 
performance. In practice, TZb is determined by directing laser pulses containing in average one photon at a given 
repetition rate on the detector and measuring the rate of counts, so in our model it is described by the term j£ . 
Analogously, TZd is calculated as the rate of counts in the absence of any input signal, so it is given by Jq D ^ ■ 

To do the comparison we need to set the values for our model free parameters: loq, 9, T = 7T and n. For simplicity 
we shall express the frequencies uiq and to in terms of respective wavelengths Ao and A. Thus we are left with two 
parameters, A and b, where b plays the role of the BP. Our general model cannot determine the BPx6 dependence 
for every kind of detector; nevertheless, one may argue that the BP and b must be proportional to each other, since 
for zero BP one should also have b — 0, because in this case the detector would be turned off. Here, we do not need 
to know the exact dependence of the BP on b provided we determine the breakdown value bs, corresponding to the 
BP breakdown at resonance, and take it as a measure of b. 

After numerical simulations we have chosen values of the free parameters that reproduce qualitatively the common 
experimental behavior [|__, ____, ___| [___[ and lie within the applicability region of the model: A = 500 nm, g = 10 11 Hz, 
T = 5 x 10 5 and n = 10~ n , so bs ~ 380, as shown in figure 1. Moreover, we verified that below bs both TZb and 
TZd have approximately linear dependence on b, in agreement with our qualitative arguments. In figure 2 we plot TZb 
for two different values of b as function of the light wavelength, where we see a good agreement with experimental 
results [___| and can check that TZb is proportional to b. We also confirmed numerically that TZd does not depend on 
the field wavelength, as expected. 

( E) _____ 

We verified that for the chosen parameters the emission terms [J„ and further terms in Eq. (|2"Tj) ] are at least 
10 orders of magnitude smaller than the dark counts term and even more for bright counts term in Eq. I|21[) . This 
confirms that detectors do not emit photons into the field. Intuitively, the emission of photons by the detector would 
be possible only at temperatures much higher than room temperature through black body radiation, which is not the 
case in experiments. Thus, in practice one is dealing only with bright and dark counts terms that act on the field 
simultaneously every time a count is registered, so the QJS takes the form 

Jp = diag (j B + Jd)p , (22) 
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FIG. 3: Normalized bright counts term as function of n in 
di-log scale at breakdown bs- At resonance (A = 500 nm) we 
have j3 ss 1/2 and far away from resonance (A = 1000 nm) 
(3 S3 [see Eq. ([24])]. 



FIG. 4: Normalized dark counts term as function of n at 
breakdown &b for different field wavelengths and the same 
graph in di-log scale in the inset. 



where diag means diagonal terms in Fock basis. In figure 3 we show the dependence of the normalized bright counts 
term J„ / J\ on n in di-log scale (for better visualization we joined the points). We note that near and far away 
from resonance - see the values of A in the caption - one has nearly polynomial dependence (linear in di-log scale) 

j( n B ) « j[ B ^n- 2 P = TZ B n"^ (23) 

with (3 R3 1/2 at resonance and (3 ss far away from it. Thus, in these cases one can write the operator dependence 
of bright counts as 

JbP = K b (n + l)~ P apa) (n+l) _/3 , (24) 

thus recovering the E- model with (3 ss 1/2 at the resonance and SD-model with (3 ss far away from it, although the 
values oUZb are different in each case (see figure [5]). 

The normalized dark counts term Jn°^ /Jq D ^ is shown in figure 4 in linear scale and in di-log scale in the inset. Out 
of resonance is almost independent on n, so in this case Jn S3 TZd = const. At the resonance, for n — one 
gets Jq D ^ = TZd and for n > we have J&l &£-K D - rT 2 ^ , where (3 ss 1/2 and £ is a number less than 1 (TZd does 
not depend on the wavelength). This means that at resonance the dark counts are suppressed in the presence of light. 
This happens because they occur when the detector in the ground state is excited by its intrinsic processes; however, 
at resonance the rate of excitations by photons from the field is much higher than by intrinsic processes, so the dark 
counts "have no time" to appear. Therefore, the operator form of the dark counts term in Eq. (|21[) is 

3 DP = TZ D [AoMo + £ArT V _/3 A] , (25) 

where Ao = |0)(0|, A = 1 — Ao and at resonance we have the E- model with (3 ss 1/2 and £ < 1. Far away from 
resonance we recover the SD-model with (3 ss and £ = 1. 

Thus we have exposed our microscopic model for the photodetector and showed that when one is concerned about 
the photodetector behavior, the model agrees with experimental data. Still, the only way to verify whether the formal 
expressions of QJSs resulting from the model hold in practice is to perform photocounting experiments and compare 
the outcomes, such as photocounts or WT distributions, to the model's prediction. In the next sections we shall treat 
this issue for realistic situations of detectors with non-unit quantum efficiency and non-zero dark counts rate. We 
also discuss possible measurements that may permit to discern between the E-model and the SD-model, even in the 
presence of low efficiency and dark counts. 
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IV. INCLUDING THE NON-IDEALITIES 



We consider a model for a photodetector with non-unit quantum efficiency (QE) and a finite dark counts rate. 
In [48| we have also considered the effects associated to cavity damping and the detector's dead-time, and showed 
that they are not crucial when compared to the QE and dark counts. Moreover, as the QJS (JTJ) is an unbounded 
superoperator, the inclusion of dead-time effect into the CPM leads to some inconsistences within SD-model, such as 
non-normalizablc photodetection distribution. On the other hand, the E-modcl is free from such problems. 



A. SD-model 



We consider a free electromagnetic monomodal field of frequency w, enclosed in an ideal cavity together with a 
photodetector. The unconditioned time evolution (UTE) of the field in the presence of the detector, i.e. the evolution 
when the detector is turned on but the outcomes of the measurements are disregarded (not registered), is described 
by the master equation 0, |24j 

p = —iu [hp — ph) — — (jip + ph — 2Ap^J , Ap = dpa\ n — a) a. (26) 

The first term on the RHS stands for the free field evolution and the second describes the effect of the detector on 
the field due to their mutual interaction. The parameter 1Z = TZb (we omit the subscript B to simplify the notation) 
is the field-detector coupling constant, roughly equal to the ideal counting rate. 

To describe photodetection with QE r\ and the dark counts rate IZd (d is the ratio between the dark counts rate 
and the ideal photon counting rate), we take (3 = and £ = 1 in Eq. (|23|l 



Jp = K (r)A + dj p. (27) 



The first term within the parenthesis describes the absorption of a photon from the field with probability per unit 



time Tr 



rfRAp 



rflin, where n is the field mean photon number - this means that the detector "sees" all the 

photons. The second term describes the dark counts, so after a detector's click the field state becomes a mixture of 
two possible outcomes: either a photon absorption or a dark count. 

The no-count state pq = Stp , where St is the no-count superoperator, obeys Eq. (|26p when one subtracts the term 
(f27|) on the RHS (see 0,111). Moreover, as we are interested in calculating probabilities, we shall disregard phase 
factors exp(±iwni), since they are canceled in the trace calculation. So the evolution equation of p s is 

TZ 

p s = -— {npg + p s h) +KvAp s -Kdp s , v = l-r). (28) 
We solved this equation in [1^], obtaining 

S tPo = e- dnt E t (e v ^ A p ), E tP = e~ nM / 2 pe- nM /\ & = 1 - e~ m . (29) 

The field UTE superoperator T t , defined as the solution to Eq. (|26|) . is naturally given by setting d = 77 = in 
Eqs. (|28|) and (f29|) . We introduced in Eq. ([29]) a compact notation for the infinite sum in terms of the exponential 
superoperator: 

exp (v4> t A)p = £ ^f-a l Po (at)' . (30) 
;=o 

The m-counts superoperator Nt(m), that describes the field state after m registered counts (whatever real or dark 
ones) in the time interval (0,£), and whose trace gives the probability for this event is 

rt rt m rh. 

N t (m)p= / dt m / dt m -i--- dtiS t -t m JS tm -t m _ 1 J • • • JS tl p, (31) 
Jo Jo Jo 

and after some manipulations [H| it reduces to 

7 . (dut+^Ar m , 

N t {m) = S t : . (32 

ml 
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A simple manner for contrasting the predictions of the model to the experimental data is by looking to the lower 
photocounts moments 



m t 



^ mTr N t (m)p 



m=0 



dlZt + nnip t 



(33) 



m(m - l) t = {dntf + 2nndllt<p t + (n<j> t ) 2 n(n - 1), 



(34) 



where n and n(n — 1) are the factorial moments of the initial density operator. 

Another measurable quantity we consider is the waiting time distribution. It describes the probability density for 
registering two consecutive clicks separated by the time interval r, provided the first one occurred at time t. Its 
non-normalized form is 



W t (r) = Tr JS T JTtp 



and the mean WT is 



p8 p8 

= Af- 1 / drW t (r) r, J\f = drW t (r) , 
Jo Jo 



(35) 



(36) 



where 9 is the time interval during which one evaluates the averaging in experiments. As to be shown in section [V] 8 
is an important parameter due to the presence of dark counts. We obtained the expression for W* (r) in [48| : 



W t (r) - e- dnr { V 2 e- K(2t+T ^Y + vde- nt (l + e^^f + d 2 ^) 



(37) 



where 



i — k 



Pn = (n\p\n). 



(38) 



B. E-model 



We now repeat the same procedures for the E-model with the QJS 

Jp = K{ife + d) Pl ep = E-pE+. (39) 

For simplicity we considered a simplified form for the dark counts term, analogous to the one we used in SD-model, 
but different from the one given by Eq. (|25|) . The probability per unit time for detecting a photon is r]lZ(l — po), 
where po = (0\p\0), so the detector "sees" whether there is a photon in the cavity. In principle, the parameter 1Z 
is different from the one in SD-model, but here it will be always clear which one we are dealing with. The field 
UTE is described by an equation similar to Eq. (|26p. obtained by doing the substitution {&,at} — > {E-,E + } in the 
non-unitary evolution [second term on the RHS of Eq. ([26]) ]. So the no-count state p s obeys the equation 

Ps = ~J + Ps A) + Kveps - a"Jlp s , (40) 

[similar to Eq. (|2"8|) ] where A = E + E- = 1 — Ao, Ao = |0)(0|. Since we are going to calculate probabilities, it is 
sufficient to write out just the diagonal form of the no-count superoperator in the Fock basis, given by [48| 

P t = e - nt ^- v ' e \ (41) 

where the dot • stands for any density operator. Once again, the functions of the superoperator e should be calculated 
as power series. 



S t 



-dm 



i 1 - P i 
Pf+Ao- * • A 

1 — ve 
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The rn-counts superoperator is 
N t (m)- = e 



-dlZt 



, * 1 A \ f, ( Jt ) m A 1 (<WW) m 1 

1- A - ■ A Pt^-i- +A - -- r— • A o 



1 — ve 

'dTZt + rjellx 



m— 1 



1 — ve m). 



•Ao 



(42) 



' 1 — ve J ' (m — 1)! 
where the last term is zero for m = 0, and the expressions for the initial factorial photocounts moments read 

m t = dKt + rin(l-E 1 ), (43) 



m(m - 1), = {dKt) + 2nndKt (1 - Si) + ff n{n - 1) (1 - Q) - 2nTZtE 



where 



S fc = ^Tr 

n 



1 -£ 



n(n — 1) 



=Tr 



Using Eq. (f35|). the WT distribution is found to be 



1-e 



P 4 °P 



, P t ° = P t (« = l). 



(44) 



(45) 



W t (r) = e- dKT <^ (ftd) 2 [l - Tr(P t °p)] + Tr[(JP T + 7edA ( 



1-Pr 

1 — 7J£ 



A )JP«p] 



(46) 



V. SOME SCHEMES FOR VERIFYING CPM 

Guided by experimental data (5lj | we chose the following numerical values for the model parameters: 77 = 0.6 for 
the QE and d = 5 • 10~ 3 for the dark counts rate (normalized by the ideal counting rate). We do not attribute 
any fixed value to 1Z since our analysis will be given in terms of the dimensionless IZt (t being the time). As many 
photodetection quantities were reported in different contexts [U, [ll], [ID, [l(| [T3, Hil HH Ell , we shall consider few of 
them that could help to decide between the SD- or the E- model. 

We first analyze the counting statistics. In figure [S] we plot ro t as function of TZt for both models for two values 
of the initial mean photon number, n — 50 and 100. Initially, rh t increases steeply due to photons absorption, and 
after some time the growth turns linear with much smaller slope due to the dark counts. We call the time interval 
during which the photons are absorbed (representing the duration of the steep increase in the number of counts) the 
effective counting time tE- In the E- model is proportional to the initial average photon number, contrary to the 
SD-model [as seen from the figure [5] and formulae (f33|) and (|43|) ]. So the experimental analysis of the dependence 
of tE on n seems to us a feasible way for verifying which model could hold in practice, because, according to the 
SD-model, 4g does not depend on n. Moreover, one could also check the validity of each model by verifying whether 
fht depends on the initial field state: in the SD-model it is independent of the field state, while in the E-model frit 
is quite sensible to it: in figure [5] one sees a notable difference between thermal and coherent states, although not 
so much between number and coherent states. This can be explained by a great difference in the values of Mandel's 
Q- factor [52j characterizing the statistics of photons in the initial state: it equals —1 and for number and coherent 
states, respectively, whereas it is very big {Qth — n) for the thermal states with big mean numbers of photons. 

Now we analyze the normalized second factorial moment 



K t = m(m — l)./m t 



(47) 



for the same initial states with mean photon number n — 50. For the number and thermal states Kt as function 
of TZt is shown in figure [SJ and for the coherent state we get K t — 1 , so it is not plotted. In the asymptotic time 
limit and for non-zero dark counts rate, the same value — > 1 holds for both models, however the transient is 
model dependent. In the SD-model without considering dark counts Kt is time- independent, Kt — n(n — l)/n 2 (n 
and n(n — 1) correspond to the initial field state), nevertheless it depends on the initial field state: Kt = 2 for the 
thermal state and K t = 1 — 1/ ft for the number state. By including the dark counts in the analysis this constant 
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400 



FIG. 5: fht in E-model for coherent, number and thermal 
states (indicated in the figure, the lower curves are labeled 
analogously) as function of time: the lower curves corre- 
spond to h — 50 and the upper - to n = 100. In the inset 
we plot fht for the SD-model, which is independent from 
field state. 



FIG. 6: K t , Eq. flU), for SD- and E- models (as indicated 
in the graph with abbreviations) for the number state (and 
the thermal state in the inset) for h — 50. For the coherent 
state one has K t — 1 at all times for both models. 



behavior is slowly modified as time goes on, sec figure [5] In the E-model in the absence of dark counts K t starts at 
the value 

liu lKt = ^J^= 1 -^, (48) 
[Tr(ep)] 2 (l- Po f 

which is exactly 1 for the number state and very close to 1 for the thermal state with the chosen values of h. With 
the course of time, Kt attains the same values as for the SD-model (for respective initial field states) when all the 
photons have been counted. By taking in account the dark counts, such a behavior is slightly modified, yet it is quite 
different from the behavior in the SD-model, as shown in the figure [6] This is another possible manner for verifying 
the applicability of SD- or E- models. 

We now turn our attention to the WT analysis. It is important to define the time interval over which we do the 
average: if one has non-zero dark counts rate, then by performing the average over a very large time interval, we shall 
always get for the mean WT the value f ~ (IZd) 1 , which is nothing but the mean time interval between consecutive 
dark counts. Since experimentally the average is done over finite time intervals, we shall proceed in the same way: the 
mean WT for initial times, when the photon number is significative, is roughly {qTZ)" 1 (because vfR. is the effective 
counting rate), so we shall take the average over a time interval 9 — 10(?77?.) _1 . This means that if one does not 
detect consecutive counts within the time such a measurement will not contribute to the average. In an ideal case 
this procedure is not necessary because the probability for registering consecutive clicks separated by a large time 
interval is zero. 

In figure [7] we plot the mean WT for the SD- and E- models, for the number and thermal initial states (for the 
coherent state we obtain a curve almost identical to the one for the number state) with n — 100 as function of the 
mean photon number in the cavity at the moment of the first click, 



Ncav = Tr 



ne nt for SD-model 
nSi for E-model. 



(For completeness, in the inset of figure [7] we plot Ncav as function of IZt for both models.) For the E-model, we 
see that when Ncav becomes less than 1, the WT starts to increase visibly due to the dominance of dark counts, 
which are much more rare events than absorption of photons. This is a drastic departure from the ideal case, in 
which no counts occur after all the photons have been absorbed, so the mean WT saturates at the inverse value of 
the counting rate, as shown in 21]. Moreover, from figure [7] one verifies that as long as there are photons in the 
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FIG. 7: Mean WT f t as function of N C av for the 
number (N) and thermal (T) states for SD- and E- 
models. While there are photons in the cavity f t is 
constant for the E-model, but increases substantially 
100 10 1 0.1 0.01 1E-3 1E-4 with time for the SD-model. In the inset we plot 

Nq^w Ncav as function of TZt (assuming the same 1Z for 

both models) for these states (in the SD-model Ncav 
is state independent). 

cavity the mean WT is nearly time-independent within the E-model (and truly independent in the ideal case [2ll|). 
and does increase substantially in time for SD-model. This is another notable qualitative difference that could be 
verified experimentally. 

VI. SUMMARY AND CONCLUSIONS 

We presented a microscopic model for a photodetector modeled as a 2-level quantum sensor plus a macroscopic 
amplification mechanism. Using the quantum trajectories approach we deduced a general QJS describing the back- 
action of the detector on the field upon a photocount and showed that it can be represented formally as an infinite 
sum of terms. In that sum we have identified the terms corresponding to the bright counts (photoabsorptions), the 
dark counts and emission events, each one occurring with its respective probability. Adjusting the free parameters 
of the model to fit experimental data, we showed that the emission terms can be disregarded in realistic situations 
since their contribution becomes insignificant, so the QJS consists effectively only of bright and dark counts terms. 
Moreover, we have simulated the experimental behavior of the counting rates and the Signal-to-Noise ratio, showing 
the breakdown phenomenon. We have also showed that with the detector operating near its breakdown bias one can 
engineer the QJS by modifying the wavelength of the field. In particular, one recovers the QJSs proposed previously 
ad hoc: at resonance one gets the E-model, and far away from it the SD-model is identified. 

We have also generalized the continuous photodetection model through a quantum treatment of non-ideal effects that 
are ubiquitous in experiments. We derived general expressions for the fundamental operations in the presence of non- 
unit quantum efficiency and dark counts, and calculated explicitly the photocounts and the waiting time probability 
distributions for initial coherent, number and thermal field states. By calculating the first and second factorial 
moments of the photocounts and the mean waiting time, we showed that in standard photodetection experiments one 
could check the applicability of the QJS of SD- or E- models. Namely, we indicated three different ways for revealing 
the actual QJS: (1) quantitatively, by studying the time dependence of the normalized second factorial photocounts 
moment. Qualitatively, we showed that the models can also be distinguished by measuring: (2) whether the effective 
detection time depends on the initial average photon number in the cavity and (3) whether the mean waiting time is 
modified as time goes on. Still, if the experimental data would depart significantly from the theoretical predictions 
one should reconsider both models and try to look for alternative mechanisms to reproduce the outcomes. 
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